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. . . Abstract 

Interactions are explored through the observation of the dynamics of particles. On the classical 
level the basic underlying assumption in that scheme is that Newton's second law holds. Relaxing 
the validity of this axiom by, e.g., allowing for higher order time derivatives in the equations of 
motion would allow for a more general structure of interactions. We derive the structure of 
interactions by means of a gauge principle and discuss the physics emerging from equations of 
motion of higher order. One main result is higher order derivatives induce a zitterbewegung. As a 
consequence the main motion resulting from the second order equation of motion is rather robust 
against modifications. The gauge principle leads to a gauge field with the property of a space 
metric. We confront this general scheme with experimental data. 
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1 Introduction 

The most basic approach to the mathematical description of nature is provided by Newton axioms 
PQ. The most important of them state that (i) there are inertial systems, (ii) the acceleration of a 
body with respect to an inertial system is given by 
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f — . . mx = F(x,x,t) , (1 

t^J- ■ where m is the (inertial) mass of the body, x its acceleration and F the force acting on the body 
which may depend on the position and the velocity of the particle, and that (iii) actio equals reactio. 

Leaving aside the fundamental and still unresolved problem of how to really define an inertial 
system (see, e.g., [2]), the second Newton axiom is a tool to explore the forces and, thus, to measure 
the fields acting on particles. The electromagnetic field, for example, can be explored and measured 
through the observation of the acceleration of charged particles under different conditions (different 
initial conditions, different charges, etc.). Thus, by means of the dynamics of the form (pQ) the 
electromagnetic interaction and, in principle, all other interactions are defined. In many cases one 
uses quantum equations of motion like the Schrodinger or Dirac equation which, via the path integral 
approach or the Ehrenfest theorem, for example, are also based on ((T|). All these considerations 
extend to relativistic equations of motion. 

From this observation it is clear that what one defines as interaction or as the corresponding 
force field depends on the basic structure of the equation of motion. If, for example, the equation of 
motion is of higher than second order in the time derivative, then interactions could be introduced 
in a different way and, thus, can have a different structure as we will show below. As a consequence, 
it is very important to have an experimental basis for the fundamental equation of motion to be of 
second order in the time derivative. 

Equations of motion of higher order are related to a different initial value problem: one needs 
more initial values beyond the initial position and initial velocity. In physical terms this means that 
with respect to the time coordinate the equation of motion is more non-local than the ordinary 
second order equations. The extreme case that (formally) an infinite number of initial values are 
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needed is related to equations of motion with memory as, e.g., generalized Langevin equations where 
the force equation possesses an additional term of the form f Q C{t — t')x{t')dt' , see, e.g., [3]. 

Also within quantum gravity scenarios it might be expected that the effective equation of point 
particles may contain small higher order time derivatives. In fact, since space-time fluctuations in 
the sense of, e.g., fluctuations of the metrical tensor, yield Langevin-like equations of motion also 
quantum gravity scenarios naturally are expected to lead to effectively higher order equations of 
motion where the additional higher derivatives probably scales with, e.g., the Planck length. 

Higher order derivatives in equations of motion occur in effective equations which take back 
reaction into account. One example for that is the Abraham-Lorentz equation for charged particles 
taking into account the electromagnetic waves radiated away [I] or the radiation damping equation 
in gravity where the emitted gravitational waves are taken into account, [S]. In the electromagnetic 
case the leading term is a third time derivative which leads to unphysical runaway solutions what still 
is an unresolved problem. However, these equations discussed in relation with radiation damping 
are no fundamental equations, they are effective equations emerging from the fact that one no longer 
regards the particles as test particles. Here we are only interested in the fundamental equations of 
motion. 

Since spin is some element of non-locality it is not astonishing that also the dynamics of spinning 
particles effectively can be described by means of a higher order theory [JjJ. 



Pais and Uhlenbeck [7] studied higher order mechanical models as a toy model for discussing 
properties of higher order field theories where higher order derivatives came in naturally in noncom- 
mutative models [8j or are introduced in order to eliminate divergences, see, e.g., {9j. Recently, it 
has been shown in [TU] that the energy of the Pais-Uhlenbeck oscillator is bounded from below, is 
unitary and is free of ghosts, see also [llj for further studies in this direction. 

In the following we will consider fundamental higher order equations of motion which can be 
derived from a variational principle. In order to be able to confront these modified equations of 
motion with experimental data one first has to investigate the structure of interactions. This will 
be done by using a gauge principle for a second order Lagrange formalism (which in principle can 
be applied to Lagrangians of all orders) . The solutions of the equations of motion coupled to these 
generalized gauge fields show that essentially only a zitterbewegung is introduced by the higher order 
derivatives showing that the standard equations of motion are rather robust against addition of 
higher order terms. We also discuss the experimental possibilities to search for effects related to 
higher order time derivatives. As an interesting by-product of the corresponding higher order gauge 
formalism we obtain the standard spatial metric as a gauge field. 

2 Lagrange formalism 

First we will use the Lagrange formalism in order to introduce interactions into a theory containing 
higher order derivatives. We introduce interactions by means of a gauge principle. We set up our 
notation by repeating shortly the standard first order formalism and then apply a gauge principle 
to a second order Lagrange formalism. 

2.1 First order formalism 

A Lagrange function of first order 



We obtain the same equation of motion from two Lagrange functions if they differ by a total time 
derivative only 





yields the Euler-Lagrange equation of motion 



= VL - — V±L. 



(3) 




(4) 
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where the function / is allowed to depend on t and x only. We can expand the total time derivative 

L'(t, x, x) = L (t, x, x) + dtf(t, x) + x- Vf(t, x) (5) 

We can now invoke the gauge principle which prescribes the replacement 

dtf{t,x)^-q<l>(t,x), Vf^qA(t,x), (6) 

where q is the coupling parameter (charge). The new function (fi(t, x) and A(t, x) are the scalar and 
vector potential of the Maxwell theory. Then the Lagrangian coupled to these potentials reads 

L'(t, x, x) = Lo(t, x, x) — qcj)(t, x) + qx ■ A(t, x) (7) 

With the choice Lo(t,x,x) = \mx 2 we obtain the standard Lorentz force equation of a charged 
particle moving in an electromagnetic field. 

Now we generalize this approach to higher order Lagrange functions. 



2.2 Second order formalism 

In the second order formalism we consider Lagrange functions of the form 

L = Lo(t, x, x, x) (8) 
from which we obtain the equation of motion 

where V a denotes the gradient with respect to the variable a. Also in this case we obtain the same 
equation of motion from another Lagrange function if it differs from the original one by a total time 
derivative of a function / only. This function, however, now may depend on the velocities x 

L(t, x, x, x) — > L'(t, x, x, x) = Lo(t, x, x, x) + x > ®) • (10) 

The expansion of the total time derivative gives 

L'(t, x, x, x) = Lo(t, x, x, x) + dtf(t, x, x) + x ■ V/(t, x, x) + x ■ Vi/(i, x, x) . (11) 

The question now is how to employ the gauge principle. If we replace, e.g., dtf(t,x,x) by a 
function 4>(t, x, x) then this function cannot describe an external field since it would depend on the 
velocity. A given external field should be given per se and should not depend on the state of motion 
of a particle. The properties of the external field should be independent of whether the particle is 
moving through it or not. 

One way to introduce functions depending on time and position only is to assume that the 
function f(t,x,x) is polynomial in the velocity. That means 

N 

f(t,x,x) = Y,h...i n (t,x)i ll ---x ln - (12) 

n=0 

In this case we can regard the functions fi- L ...i n {pi x ) = f(i 1 ...i n )^Pi x ) as § au S e functions of an exter- 
nally given interaction. For this setting we obtain for the new Lagrange function yielding the same 
equations of motion 

N 



L'(t,x,x,x) = L (t,x,x,x) + y^d t fi 1 ...j n (t,x)x 11 ■ ■ ■ x ln 

71=0 

N N 



n=0 n=l 
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The gauge principle now allows to replace these gauge functions by the gauge fields 

dtfh...i n (t,x) -> -q n 4>i 1 ... in (t,x) 

d i fi 1 ...i n (t,x) -► q n A ih ... in (t,x) (14) 
nf h ... in {t,x) -> q n ipi 1 ...i n (t,x), (15) 

where the g n are the coupling parameters to these n th rank potentials. The symmetries of these 
gauge fields are 

&i...i„(t, = 0(u...j„)(t, as) , ^...^(i, x) = A (iil ... in) (t, x) , Ai...i n {t, x) = ip {h „ An) (t, x) , (16) 

that is, all gauge fields are totally symmetric. These fields transform under the generalized gauge 
transformations as 

qn<f>h...i n (t,x) -> q n <l>i 1 ...i n (t,x) = q n <f>h...i n (t,x) - d t f il ... in {t,x) 
q„A ih _.. in (t,x) -»• q n A' iilin (t,x) = q n A iil ... in (t,x) + d( i f il .. An) (t,x) (17) 
qn^h...i n (t,x) -> q n i)' h _ in (t,x) = q n ij) h _ in (t,x) + nf il ... in (t,x) . 

For n = they reduce to the ordinary Maxwell gauge transformations. 

As a consequence we obtain the Lagrangian of a particle coupled to the new interaction gauge 
fields 

N 



L'(t,x,x,x) = L (t,x,x,x) - ^2q n (f)i 1 ...i n (t,x)x l 

n=0 

N N 
+ Y qnA ih ...i n (t, x)x l x h • • • x in + Y q n i>h...i n {t, x)x h x i2 • • • x in . (18) 



n=0 n=0 

The equations of motion read 

. d dV d 2 dV 
= djL'-— — - + 



dt d& dt 2 dxi 
d dL d 2 dL c 
Jt~d& + dt 2 ~dx~i 

+ Y '/»^- , » •••••• (*> x)x i x il ■ ■ ■ x in + Y QnOjAi...i n {t, x)x h x i2 ■■■x i 



d dL d 2 dL A < . in 

" w77to + w72 Tto " 'i^hor...,,, (t, *)* 1 • • • * " 



n=0 

N N 



n=0 n=l 

N N 



j t ( - Y 'I ""° J >.>.-■:- (*, x )x i2 ■ ■ ■ x in + Y ( n + \'<ln-y,, (t, x)x h ■■■x i 

' n=0 n=0 



JV \ 
+ J^(n - l)qn^i lji3 ...i n (t, x)x ll x t2 x 13 ■ ■ ■ x ln 

n=l J 

+ ^2 [Y^ j i,..i n (t,x)x^...x i A . (19) 



d 2 /N 



The additional gauge interaction can add a time derivative of at most third order. The principal 
part of the differential equation remains unaffected. 

Below we typically will use the most simple second order Lagrangian without interaction 

L (t,x,x,x) = ^x 2 - ^x 2 . (20) 

The parameter e has the dimension kg s 2 . If one assumes that this additional term has emerged from 
influences of quantum gravity then it might be natural to identify it with e ~ MpjTpj ~ 10~ 95 kg s 2 
which is extremely small. 
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5S = I [VL-—V ± L + —;VxL)5xdt 



3 Noether theorem 

Also for higher order Lagrangians conservation laws can be obtained from the variational principle 
if we allow non-vanishing variations at the initial and final points. The variations are as usual 

t = t + r(t) (21) 
x(t) = x(t) + Ax(t) . (22) 

One should bear in mind that here the x need not to be variables of the configuration space. For 
these general variations, the variation of the action is 

SS = S-S = L(x(i),x(i),x(i),t)dt- f ' L(x(t),x(t),x(t),t)dt. (23) 
Jti Jti 

Proceeding as in the first order case described in textbooks we arrive at 
rta d d 2 

V & L ■ j t (Ax - tx) + (v ± L - j^xL^j ■ (Ax - tx) +rij ' . (24) 

If the equations of motion are fulfilled, and if the action is invariant under the variations (|2 1I22|) . 
then we obtain the conserved quantity 

p 2 ■ — (Ax — tx) + (p 1 — p 2 ) ■ (Ax — tx) + tL = const , (25) 
at 

where we defined the momenta 

p 1 = V±L , and p 2 = V £ L . (26) 

If the action does not vanish but, instead, changes with a total time derivative of a function 
F(x(t),x(t),t), then the equations of motion in terms of the Euler-Lagrange equations do not 
change and we have a modified conserved quantity 

p 2 • -j- (Ax — tx) + (p 1 — p 2 ) ■ (Ax — tx) + tL + F = const. (27) 
From this general Noether theorem we derive the following conserved quantities: 

3.1 Momentum conservation 

At first we consider the transformations r(t) = and Ax(t) = const. We obtain the conserved 
momentum 

P = Pi — p 2 = const. (28) 

This may also directly be inferred from the Euler-Lagrange equations of motion ([9]). For the 
Lagrangian (f2"UJ) we get 

P = mx + ex = const. (29) 

3.2 Energy conservation 

Next we consider the transformations r(i) = To and Ax(t) = 0. The corresponding conserved energy 
is 

E = p 2 ■ x + (pi — p 2 ) ■ x — L = p 2 x + P x — L = const . (30) 
For the Lagrangian (|2U|) we obtain 

E = -mx 2 + -e (2x ■ x - x 2 ) . (31) 
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3.3 Angular momentum conservation 

In the next example we assume that the action is invariant under the transformations r(t) =0 and 
Ax = A<p x x. This corresponds to a rotation and the corresponding conserved angular momentum 
is 

L = x x (p 1 — p 2 ) + xxp 2 = xxP + xxp 2 = const. (32) 
For the Lagrangian (|20f) we obtain the conserved angular momentum 

L = x x (mx) + e(x x x — xxx) = xxP — ex x x = const. (33) 



3.4 Proper Galileo transformation 

At last we consider the Galilei transformations r(t) = and Ax = Avt for the Lagrangian (|20p . 

where Av is assumed to be very small. The term |mi 2 changes by a total differential of mx ■ Av 
so that we obtain the conserved quantity 

C = p 2 + (Pi ~ P 2 ) t ~ mx ( 34 ) 
from which we deduce a uniform motion 

P2 1 Pi ~ P2 . , /orN 

x = 1 1 + xq . (35) 

m m 



4 The most simple gauge model, n = 
4.1 Equation of motion 

The most simple case with non-trivial dynamics is given for the special case that / in (|lip is a 
function of t and x only. Then we have 

L'(t, x, x, x) = Lo(t, x, x, x) — qcj)(t, x) + qx l Ai(t, x) . (36) 

For Lq from (|2(J[) the equations of motion read 

e x +mx = qE(t, x) + qx x B(t, x) , (37) 

where E and B are the electric and magnetic field derived as usual from the scalar and vector 
potentials 4> an d A. This equation of motion is the standard one with a small additional forth order 
term. 

We may simplify even further by taking a vanishing magnetic field, B = and a constant electric 
field E(x) = Eq = const, 

e x +mx = qEo . (38) 
This is the equation we like to solve now. 



4.2 Particle motion 

The first two time integrations are easily performed and give 

ex + mx = ^E {t - t ) 2 + a(t - t ) + b , (39) 

where a and b are integration constants. Next we introduce a new function x through x = xq + ex 
with xq = ^Eo^t — to) 2 + -^a{t — to) + ^b which represents the solution of the corresponding 
equation of motion without the forth order term. Since e will be assumed to be small, the deviation 
from the standard solution should be small, too. Therefore we introduced an extra e in front of x. 
The differential equation for x then reads 

■^ + —x = -S-Eo. (40) 
e me 
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A further substitution x = x + ^^Eq yields the equation for a particle in a harmonic potential 

x + —x = 0, (41) 
which, according to the sign of e, possess the solution 

x = A cos (ut) + B sin {ut) for e > (42) 

x = A\ cosh (ut) + A 2 sinh (art) for e < (43) 



for some amplitudes A, B, A\, and A 2 and with uj = \Jmje. 
As a consequence we arrived for e > at the solution 

x(t) = —E (t - t ) 2 + — ait - t ) + — b + eA cos (ut) + eB sin (ut) - ——E . (44) 
2m m m mm 

In the limit e — ► we obtain the standard solution. The velocity is 

x(t) = —Eo(t — to) H a — ^JmeA sin (ut) + ^frneB cos (ut) , (45) 

mm 

which also approaches the standard expression for e — » 0. The acceleration 

£c(£) = — £o — m A cos (wt) — mB sin (ut) , (46) 
m 

however, has a large fluctuating term of order 1. 

For very small positive e the additional term in the path (j44f) is very small but fast oscillat- 
ing zitterbewegung. One may turn the above result into a positive statement: At least within a 
Lagrangian approach and assuming that higher order terms are small, the paths originating from 
a corresponding higher order modification are rather inert against these modifications. The mean 
orbits behave like orbits given by second order equations of motion. We will present some ways to 
experimentally search for these fundamental oscillations in Section [6j 

This result may be extended to the case of slowly varying arbitrary electromagnetic fields. 
The equation of motion for an arbitrary electromagnetic field (|37H may be attacked through the 
substitution x = ex + xq where xq is assumed to solve the equation of motion without the forth 
order term. If we assume that the force F(x) = qE(x) + qv x B(x) is very smooth and that the 
deviation ex is very small (that is, if x ■ V.F -C mx and can be neglected), then we obtain 

xq +e x +mx = . (47) 

This can be integrated twice 

xq + ex + mx = at + b , (48) 
where a and b are two integration constants. Inserting the equation for xq yields 

Tfl \ ' I 

x-\ x = F(x ) + -at + -b. (49) 

e me e e 

With a new variable x = x — —at + —b KF(xn) we have 

Tfl 

x + —x = 0. (50) 

Then again, for positive e, x is a fast oscillating term which adds to the standard solution. The 
total solution then is 

x(t) = x (t) + e (x(t) + —at -—b+ —F(x (t))) . (51) 
\ mm m z / 

This solution consists of the standard solution Xo(t) which is the main motion, a small position- 
dependent displacement and a small linearly growing term, and a small fast oscillating term, a kind 
of zitterbewegung. From ordinary observations, a and b should be very small. Neglecting these 
particular solutions, then the standard solution of the standard second order equation of motion 
seems to be rather robust against adding a higher order term. 
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4.3 Conserved energy 

The conserved energy in this most simple model reads 

E = y x 2 + | (2x ■ x - x 2 ) + q<t> , (52) 
which also can be obtained by multiplying the equation of motion by x and partial integration. 

5 The next gauge model, n = 1 
5.1 The equation of motion 

Here we now like to explore the model where / is a function which is polynomial of first order in 
the velocities, that is, f(t, x, x) = f^°\t, x) +f- (t, x)x % . Then the corresponding gauged Lagrange 
function reads 

L(t, x, x, x) = L {t, x, x,x) - q (f) + qox l Ai - qi^fa + q\x l x? + q\x % i>i , (53) 

where all functions depend on t and x. Here qo is the coupling related to f(°>, and qi the coupling 
related to /j. Beside the usual scalar and vector potential, <j) and Ai we have in addition two vector 
potentials fa and tpi and a tensorial potential A^j. The gauge transformations are given by 



q Q( f> - dtf® , qoK = qoAi + fyf^ 

' Vi) 



qi4>'i = <li4> - d t fi l) , qM = qiA + fi l) , qiA'^ = q x Aij + d {i f^ 



(54) 



The corresponding Euler-Lagrange equation is 



dL d dL d 2 dL 



° = ~dx^ ~ ~dt ~d¥ + dfi ~d¥ ^ 
-q di<j) - qod t Ai + q\d t fa + q\d\fa (56) 

+x 3 (diiqoAj - qifa) - (djq Ai - qidjfa) + 2q 1 d t d [j i(j i] ) (57) 
-2« + + a,A ik - a,A Jt) ^) - 1, A A^ . (58, 

where we defined Ay := A^j — d^ipjy It should be noted that the gauge terms do not produce any 
x '-terms. 

The first line (|55p is the (still unspecified) equation of motion of the free particle, the second 
line (j56|) describes a force due to two electric fields E^ := — V</> — dtA and E^ := dt4> + dfij), 
the third line (|57[) is a Lorentz-like force with the two magnetic fields B^ := V x A and := 
V x (— <f> + <9t ■?/>), which are all gauge invariant. The forth line (j58[) resembles the form of covariant 
derivative with a connection based on a second rank tensor Aij. 

Taking again Lq of the form (|2Up yields an equation of motion of forth order of the form 

e x j Sij + m gij (D ± xy - md t g l3 x^ = q Ef ] + qi E^ + (x X (q B^ + 31^% . (59) 

where we introduced an effective 3-metric 

9ij = S i:j + 2^ Aij (60) 
m 

and where the covariant derivative D± is formulated with a Christoffel symbol based on g^j. This 
metric is invariant under the gauge transformations (|54|). The second and third term on the left 
hand side can be regarded as equation of motion arising from the variation of J gij(t, x)x % x J dt. If in 
(|59[) we let e — ► then only the forth order derivative term will vanish, the other terms, in particular 
the metric, remain. 
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One should note that by means of our second order gauge principle we were able to establish a 
space-metric as an ordinary gauge field. Together with this metric gauge field we also introduced 
an interaction with an electromagnetic gauge field. Therefore, by means of one gauge procedure we 
introduced the the interaction with a metric as well as with an electromagnetic field. 

By omitting the n = part of this gauge formalism, or equivalently, by setting qo = 0, we have a 
combined gauge formalism for the metric and an electromagnetic field. We also can choose .E (1) = 
and B^ = while keeping Aij ^ 0, E^ ^ 0, and B^ ^ 0. This latter choice introduces the 
electromagnetic field and the space metric independently. 

The structure of the above equation of motion (foTJj) is 

e x +mx = K(t, x, x) , (61) 

where K(t,x,x) is a polynomial of order two in the velocities what slightly generalizes (I37D . As a 
consequence, for ordinary situations (smooth forces) and small e we again obtain some zitterbewegung 
resulting from the inclusion of higher order derivatives. Again, the main (mean) motion is described 
by the second order part of the equation of motion. 

5.2 The conserved energy 

We assume that all fields do not depend on time and determine the conserved energy. Then again 
we obtain an expression where the space metric appears at the right place. Applying ([30]) yields 

1 1 

E = -mgijx'x 3 - q <p + -e (x 2 - 2x ■ x) . (62) 

As expected, the effective metric (|6Up enters the kinetic energy. 

5.3 Relativistic formulation 

It is clear from ^ipi = —% i 4f'<Pi = —x % dt^% — x^x^djipi which holds modulo a total time derivative, 
that 

- q\x % <j>i + qix l x j Aij + q\x % in = -q\x l ((pi + d t ^i) + q\x l x j (A^ - djipi) , (63) 

what explains the definitions of and Furthermore, with x a = (t,x l ), a = 0, ...,3, and 

the 3+1 splitting 

q A a x a + qi A ab x a x b = q A + qi A 00 + (q Ai + 2 qi A i0 ) x l + q^xW , (64) 

we can reproduce, with redefinitions, the gauge ansatz (f53l) above. Therefore, the gauge field part 
in the Lagrangian (|53|) can be written in 4-covariant form 

L(t, x, x, x) = Lo(t, x, x, x) — qoA a x a + q\ A a b± a x h . (65) 

It is the last term quadratic in the velocities which, when interpreted as gauge field, requires the 
second order Lagrangian. 

6 Comparison with Experiments 

We like to describe possible experiments which may be sensitive to the higher order modifications in 
the equations of motion. We take situations where a charged particle is placed in a constant electric 
field, e.g., within a condenser. In standard theory the charged particle will be accelerated and the 
final velocity depends on the voltage only, not on the spacing between the plates of the condenser. 
To make the situation as simple as possible we consider a one-dimensional problem and take as 
initial conditions that the particle is at absolute rest at to = 0. 
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We take our general solution (|44j) (in one dimension and with to = 0) 



x(t) = —E t 2 + —at + —b + eAcos(cot)+eBsm(cot)-——E (66) 
2m mm mm 

x(t) = —Eot-\ a — me sin (cot) + B\/me cos (cot) (67) 

m m 

x(t) = —Eq — Am cos (wi) — Smsin (cot) (68) 
m 

3 3 

x(t) = eA (J^j 2 sin (tot) - eB (J^j 2 cos (cot) (69) 

and determine the parameters in terms of the initial conditions x(0) = 0, x(0) = 0, x(0) = 0, and 
x'(0) = 0. We obtain B = 0, a = 0, A = -\Eq, and 6 = 0. The solution then reads 

x(t) = ±E (It 2 + - (cos (wt) - 1)) (70) 
m \z m / 

and we also have 

x(t) = — E ( t- J— em(ut) ) (71) 
m \ V m / 

= (1 - cos (cot)) (72) 

m 

C / Tfl 

x(t) = —E <—sm(ujt). (73) 
m \ e 

Now we discuss three different measurements which outcome depend on e: (i) the time of flight 
in an accelerator, (ii) position and velocity variance, (iii) noise in electronic circuits. We always 
calculate first order effects. We also do not take into account possible effects of radiation reaction. 
Due to the modified fundamental equation of motion this needs a new derivation of this effect which 
will be planned for the future. 

6.1 Determination of time of flight 

In the standard theory (e = in the Lagrangian) we have the solution x(t) = -^Eot 2 . For a 
charged particle being accelerated through its motion through an accelerator of length L we have 
L = z^^-t 2 , where Acj) is the potential difference the particle has to traverse. Then t 2 = -^^L 2 
so that t is proportional to the spacing L and the final velocity does not depend on the distance L. 
This will be different for our higher order theory. 

We determine the time the particle needs to traverse the length L. For that we first have to 
calculate t from L = x(t): 

L = ^-Eo (\t 2 + - (cos (cot) - 1)) , (74) 
m \ 2 m J 



what is a transcendental equation. For e = we have t = to = y^g^- Therefore we make the 
approximative ansatz t = to + ti, where t\ <C to and t\ is of the order e. Solving for t\ gives the 
time of flight to first order correction in e so that 



f cos W^j - l ) (75) 



We use this time in order to calculate the velocity at x(t) = L and obtain 
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where we substituted Eq = — ^ and also inserted the velocity of the standard theory x$ = J — 2q ^f ■ 
In the standard theory x(L) does not depend on L. Here, by varying L we get first oscillations in 

the velocity due to the sin and cos terms, but also an offset 

Since e is assumed to be small, then a change in L will result in fast oscillations which probably 
cannot be resolved. Therefore, averaging over a small L interval yields 

Wi »= i »( 1+ 4^l)- (77 > 

Therefore in the mean the velocity after traversing the condenser is a bit larger. 
Rewriting the above result as relative velocity deviation 

(x(L)} - xp = e x\ 
xq Am L 2 

it is clear that one obtains good estimates for e for large velocities xq, short L and small m. Taking, 
e.g., an ion of mass of m = 60 u, a final energy of 10 MeV, a traversed distance of L = 1 m, and an 
accuracy to measure the relative velocity of 1%, then we arrive at an estimate e < 10 -40 kgs 2 . 



6.2 Interferometry 

Acceleration can be measured directly , e.g with atomic interferometry. This has been proposed first 
by Borde [12J and the today's best performance gives an uncertainty of the measured acceleration 
of Ax ~ 10 -8 m/s 2 |13j . However, while this accuracy is valid for a constant acceleration, in our 
case we have a fast varying acceleration. 
We use the phase shift 

A0 = A(uj)kxT 2 , (79) 

with our acceleration (I72p . Since that part of our acceleration we are interested in and which we like 
to detect this way is fluctuating, we have to amend the standard phase shift for a dc acceleration, 
Acp = kgT 2 , by a transfer function A(u) which has been determined in [14J. We use as charged 
particle ionized Helium and take typical values for the laser wavelength A = 780 nm, a pulse spacing 
time of T = 100 ms, an electric field strength of Eq = 10 10 V/m. An experimentally reachable 
accuracy of the phase measurement is Acj) = 10 -3 rad. 

We are interested in the largest u which we are able to detect. With the specifications given 
we can determine that u> from the condition A(u>) = 10 -25 . This gives u = 10 12 Hz which is the 
maximum frequency which effect on the phase shift we are able to detect. If we assume that nothing 
is detected this gives an estimate of 

e < H = 10~ 50 kgs 2 . (80) 

t, 



6.3 Electronic noise 

Another situation where a kind of zitterbewegung of a charged particle may induce an effect is 
electronics. As a most simple model we may assume that the particle under consideration is located 
within a capacitor. An oscillation of this particle will induce an electronic noise in the electric 
circuit. This electronic noise can be estimated by 

C(U 2 ) = m(x 2 ) (81) 

where U is the voltage. Using only the oscillating terms in the velocity we obtain 

C(U 2 ) = \e(^EX (82) 
2 \m J 
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Therefore, a modified dynamics will induce an electronic noise - beside the noise like Nyquist noise 
or a shot noise. However, due to its different characteristics, it might be disentangled from the other 
noise sources. Our noise should be a fundamental noise not depending on temperature, the finite 
number of charged particles in the electric circuit, etc. A good cryogenic noise limit is of the order 
1 nV/VTTz |15| in a wide frequency range, that is 1 nV for a measurement of duration 1 s. Taking 
as granted that no fundamental noise of this kind has been seen under the conditions of a molecular 
vacuum and cryogenic temperature, we may get an first estimate e < 1CP 68 kgs 2 , where we took 
a capacitance of 0.48 pF, a distance between the capacitor plates of 15 fim, a voltage of 1000 V. 
Taking into account a bandwith of 1 GHz, we get a more realistic estimate e < 10~ 50 kgs 2 . 

7 Conclusion and outlook 

We studied the physics resulting from a hypothetically given higher order equation of motion sub- 
stituting standard Newtons law. In order to confront such a model with experiments we first have 
to introduce the coupling to external force fields. This we employed via the standard gauge princi- 
ple. We therefore started from a higher order Lagrangian and first introduced the interaction with 
external fields through a standard gauge principle. Since it is applied to higher order Lagrangians, 
the resulting gauge field have a richer structure than in the ordinary first order formalism. Beside 
the ordinary gauge fields usually obtained in first order formalism, we also obtained a space metric 
as a gauge field. As a by-product we obtained as a second order gauge field the standard space-time 
metric. 

Then we discussed physical consequences of equations of motion with higher order time deriva- 
tives coupled to these gauge field. We solved the equation of motion for the simplest case and 
deduced observational consequences in a few very simple experimental situations. Small higher 
order terms may influence the time of flight of accelerated particles, it yields a large acceleration 
fluctuation accessible in atomic interferometry, and also a fundamental noise induced in electronic 
devices. No deviation from, Newton's second law has been observed. Very rough and preliminary 
estimates for the parameter characterizing the higher order term could be derived. However, a 
comparison with the a Planck scale version of this parameter shows that the experimental estimates 
are far away from testing quantum gravity. 

The next step is to implement a higher order theory for point particles in a relativistic context. A 
possible Lagrange function to start with might be L = m\j r) a t,x a x b + erfabX a x b , where the dot is the 
derivative with respect to some parameter along the path. We also like to implement our principle 
in higher order field theory, e.g., for scalar field Lagrange densities of the type C = r] ab d a (f>* db(f> + 
li , rf h 'rf d 'd a d c (^*d} 1 d^. As further step would be to set up some field equations for the new gauge 
fields. A first guess might be to have the usual Maxwell equations for the A a and the Einstein field 
equations for the A^. 
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